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When smooth, zero-on-average, periodic magnetic and electric fields are applied to a carbon 
mono-layer (graphene), a gap between the valence and conduction band is introduced. Here this 
gapped state is studied analytically. It is found that it does not correspond to a band insulator: 
a constant electric field induces a quantized Hall current even though the magnetic flux through 
the sample is zero and there are no Landau levels. The phenomenon is of the same type as that 
discovered by Haldane for a graphene sample in a periodic magnetic field that is not smooth, i.e. 
varies rapidly on the scale of the graphene lattice constant. The effect can be explained in terms 
of the topological theory of Thouless, Kohmoto, Nightingale and den Nijs. For the system studied 
in this paper, an explanation in terms of simple physical principles is also presented. Thus some of 
the mystery is taken out of the apparently strange phenomenon of a Hall effect without magnetic 
flux. Furthermore, Haldane's model requires control over external magnetic fields on length scales 
less than an angstrom and is therefore hard to realize experimentally. For the model studied here, 
control over external fields on length scales that are larger by two orders of magnitude or more is 
sufficient. The model is therefore more amenable to experimental realization. 



PACS numbers: 72.15.-v, 73.43.-f 

The filled conduction band and the empty valence band 
of a carbon mono-layer (graphene) touch at two incquiva- 
lent points in the Brillouin zoneji Regions of the Brillouin 
zone in the vicinity of these points are called valleys. In 
each of the two valleys, low energy excitations are de- 
scribed by a two-dimensional massless Dirac equation. 
As a consequence, electrons in graphene can propagate 
through electrostatic potential barriers. In the jargon of 
relativistic quantum mechanics, this is known as Klein 
tunnelingi^i^ When a gap is induced between the valence 
and conduction band, Klein tunneling is suppressed and 
the electronic properties of graphene change radically. 
Several distinct situations have been uncovered: 

Recently a method for producing a band insulat- 
ing state was demonstrated experimentally.'^ It involves 
growing the graphene sample on an appropriate substrate 
that induces a potential energy difference between the 
two triangular sublattices that constitute the sample's 
honeycomb lattice. In a seminal paper by Haldane,^ 
it was demonstrated that graphene's honeycomb lattice 
also supports a fundamentally different kind of gapped 
state which I will call anomalous. A signature of the 
anomalous state is a Hall effect without Landau levels 
and in a zero average magnetic field. The resulting Hall 
conductivity is a topological invariant 

In Haldane's model, the anomalous state is realized by 
the combined effect of next nearest neighbor hopping to- 
gether with a magnetic field that has the same periodicity 
as the honeycomb lattice and zero flux through a lattice 
unit cell. Because of the very short length scales (~ lA) 
on which this magnetic field must be controlled, the ex- 
perimental realization of Haldane's model does not seem 
feasible. Recently, Kane and Melet^ pointed out that the 
spin orbit interaction can generate an anomalous gapped 
state in graphene without the need for any external mag- 
netic field. In this case, since time-reversal symmetry is 
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not broken, no Hall effect with electric charge current is 
produced. Rather a spin-Hall effect occurs. However, the 
effect may be hard to observe because of the smallness of 
the spin orbit interaction in graphene. 

A gapped state can also occur when periodic mag- 
netic and electric fields that are smooth on the scale 
of the graphene lattice constant are applied jointly to 
a graphene sample. Since these fields have a periodicity 
3> lA, experimental realization is more feasible than in 
the case of Haldane's model. In this paper I show an- 
alytically that the resulting gapped state is of the same 
anomalous type as in Haldane's model. Apart from being 
experimentally relevant, the system studied here has an- 
other interesting feature. Its Hall effect can be explained 
without recourse to the topological theory. Because of 
this it is possible in the present instance to identify a 
physical mechanism responsible for the Hall effect. This 
hopefully takes some of the mystery out of the general 
but abstract topological theory,^iSiiaiiiii^ii^iiiii^ 

A gapped graphene state that results from smoothly 
varying fields was first revealed in the numerical study 
of Ref. [IB (see Fi gure 2, last frame, of that work) and 
is also mentioned in Ref. The effective magnetic 

field studied in these references is produced by lattice 
deformations. The single- valley physics of such a field 
is identical to that of a real magnetic field. The ana- 
lytical results for a single valley that is presented here 
confirm and explain numerical results obtained in Refs. 
in and |T3. A caveat is warranted though: Despite the 
close similarity between a graphene system with a real 
magnetic field and one with a deformed lattice, there is 
a fundamental difference. The gapped state produced by 
deformations is not anomalous. The reason is that for 
deformation induced fields, the anomalous behavior of 
the two valleys cancel each other, while for real magnetic 
fields, the behavior is the same in both valleys and there 
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is no cancellation. 

Several other publications deal with systems similar 
to that studied here. In Refs. [l^ and [l^ the two- 
dimensional non-relativistic Schrodinger equation with 
a periodic magnetic field that is zero on average was 
studied. It is interesting to compare this system to the 
graphene system in the limit of strong magnetic fields, 
when the magnetic length is much smaller than the pe- 
riodicity of the fields. In this regime the anomalous 
Hall effect^° was found to disappear in the case of the 
Schrodinger equation. It remains present in graphene 
as is shown here. Another relevant study^i dealt with 
graphene in a periodic electric field alone. It was found 
that a strong electric field can induce additional touch- 
ing points between graphene's valence and conduction 
bands. In the present work, only the case of weak elec- 
tric fields will be considered, so that additional touching 
points do not arise. What happens when both the elec- 
tric and magnetic fields are strong, or when the electric 
field is strong but the magnetic field weak, is a topic for 
further investigation. 

There are also at least two relevant works in the carbon 
nanotube literature. In Ref. 22 the low energy spectrum 
of a carbon nanotube in a magnetic field that is trans- 
verse to the axial direction is calculated analytically. In 
Ref. 23 the same system is considered and shown to 
support a quantized Hall current in the axial direction. 
The nanotube system can be considered a special one- 
dimensional case of the general two dimensional system 
considered below. In order to induce a gap in the nan- 
otube spectrum, an electric field must be applied parallel 
to the magnetic field. This was not considered in the 
cited references. Consequently, the Hall effect found in 
Ref. m is fundamentally different from the Hall effect 
considered here. This will be discussed in more detail 
below. However some results that hold for the system 
discussed in this work should also hold for the nanotube 
system. Where there should be agreement between re- 
sults in this paper and in the nanotube studies, agree- 
ment is found. 

The plan of the paper is as follows: In Sec. Uthe sys- 
tem to be studied is defined mathematically. In Sec. |TT] 
its zero energy eigenstates are found when the electric 
field is zero, in which case there is no gap. In Sec. Illll the 
low energy description of the gapped state is developed. 
In Sec. IIVI the topological theory of the anomalous Hall 
effect is briefiy reviewed and then applied to the graphene 
system in order to calculate the quantized Hall conduc- 
tivity. In Sec. |V]the anomalous Hall effect is considered 
again, this time without invoking the topological theory. 
In Sec. IVII an example is presented in order to illustrate 
the general results of the previous sections. Conclusions 
are presented in Sec. IVIII 



I. STATEMENT OF THE PROBLEM 

The inter-atomic distance between carbon atoms in 
graphene is approximately 1.42 A. I consider a graphene 
sample in the presence of static electromagnetic fields U 
and A that vary smoothly on this length scale. A long- 
wavelength description is therefore appropriate. This in- 
volves the Dirac Hamiltonian 

H = va ■[-iTidr + eA{r)] + U{r). (1.1) 

where <r = [axjO^y) are standard Pauli matrices, r = 
{x,y), and dr = {dx,dy). Close to the charge neutral- 
ity point, electrons in graphene obey the Dirac equation 
= H'^ . In principle two species of Dirac fermion 
should be distinguished because the low-energy spectrum 
of graphene consists of two sets of Dirac cones. How- 
ever, since the fields A and U are smooth, no scattering 
between species is possible. Furthermore, by employing 
here the valley-isotropic representation,'^ the Hamiltoni- 
ans for the two valleys are rendered identical, not only 
in form, but also in the actual values assumed by the 
parameters. An index that distinguishes between valleys 
is therefore omitted. Results obtained will apply equally 
to electrons in both valleys unless otherwise indicated. 

I consider a magnetic field Bz — dr y- A and a scalar 
potential U that are periodic in space with the same pe- 
riodicity and that average to zero. The periodicity of the 
fields define a Bravais lattice with a primitive unit cell 
denoted UC and basis vectors a^, k = 1, 2, such that 
B{r + ak) — B{r) and similarly for U . 

In what follows, any function that has this property 
will be called J7C-periodic. Furthermore, whenever the 
term 'lattice' is used below, it will refer to the magnetic 
field lattice, not the crystal lattice of the carbon atoms 
in graphene. The two lattices are entirely distinct and 
only the former is relevant for present purposes. As a 
consequence of the smoothness of A, UC is much larger 
than the unit cell of the graphene crystal lattice. Because 
B averages to zero, the magnetic flux through UC is zero. 

The graphene sample is assumed to be much larger 
than UC and only bulk effects will be considered. It is 
therefore convenient to impose periodic boundary condi- 
tions 

*(r) = *(r + r!afe). A: = 1,2, (1.2) 

where ^ 1 is an integer. All functions with this prop- 
erty will be called sample-periodic. Note that because Q. 
is an integer, any C/C-periodic function is also sample- 
periodic. 

A reciprocal lattice is associated with the real space 
lattice. The basis vectors b^-, j = 1, 2 of the reciprocal 
lattice are defined through the equations bj -ak = 2i:5j^k- 
The Brillouin zone is a primitive unit cell of the reciprocal 
lattice, and will be denoted BZ. It is convenient to define 
the Fourier components 

B^^^r^ I dr2e-^^- ''B(r), (1.3) 
Juc 
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of the magnetic field. Here fcmn = nT-bi + nb2 with m 
and n integers and P = \ai x is the area of UC. Note 
that Bqo = due to the zero-flux condition imposed on 
B and -Bm.n — B*_^ _„ because B{r) is real. 

II. ZERO ENERGY MODES 

The task is now to analyze the low-energy (i.e. 
small |£|) spectrum of H. 

For U — Q, zero energy eigenstates 4* — ((/>+, 0_) of H 
solve 

[da: ± idy + i{aa ± iay)](p± = 0, (2.1) 

with a — eA/h. Following Jackiw,^*^ I work in the 
Coulomb gauge where dr -A^ 0. Then Eq. (|2.ip . subject 
to the sample-periodic boundary conditions of Eq. (|1.2p , 
is solved by 

^±(r) =c±exp[±F(r)], (2.2) 

where c± are arbitrary constants. The real function F is 
?7C-periodic and satisfies 

dlF^eB/h. (2.3) 

In terms of the Fourier components -Bmn of the magnetic 
field, F is given by 

^(-) In^^"-^- (2-4) 

One demonstrates uniqueness of the solutions in Eq. 
(12. 2p as follows. Consider the equation for Suppose 
that 4)-{r) solves Eq. (|2.ip and is sample-periodic as re- 
quired by the boundary conditions. Define a function 
g{r) = (^(r)e^^''^. Since F is f/C-periodic, g is sample- 
periodic. The function g satisfies the Cauchy-Riemann 
equation {dx — idy)g = 0. Thus g is an analytical func- 
tion of the complex variable z = x + iy. As a consequence 
g is either constant or unbounded in some directions at 
large r. Unboundedness is incompatible with g being 
sample-periodic and therefore g must be constant so that 
4>-{r) = ge^^'^^^ is proportional to <?!>-(r) ~ c_e^^(''^ 
The uniqueness of is demonstrated in the same man- 
ner. 

The Dirac equation H'i' = E'i> with U — therefore 
has exactly two zero energy solutions which I will take to 
be 

^ , , N+ f e^^l-) \ r ^ , f ^ \ 

*+,o(r) = — 1^ j' = IT U-"<^V ' 

(2.5) 

with N± = [Jjj^ dr"^ exp ±2F(r)] 2 

ensuring normal- 
ization of ^'±,0 to unity over the sample. The two eigen- 
states have opposite sublattice polarization: In valley K, 
is A sublattice polarized while ^'_^o is i?-sublattice 



polarized. Since the valley isotropic representation is em- 
ployed, the roles of the sublattices are reversed in the 
K' valley. Hence, the state ^+,0 in the K' valley is B 
polarized while the K' state VP-^o is A polarized. The 
result is consistent with the Atiyah Singer index theo- 
rem. It states that in a given valley the difference be- 
tween the number of zero energy eigenstates that are A 
polarized and B polarized equals the total flux through 
the samplci^ Note also that because F is C/C-periodic, 
the same is true for VPi^Oi even though only sample- 
periodicity is imposed on \E'±,o by the boundary condition 
[Eq. (1121)]. 

The solutions in Eq. (j2.5p were found without invoking 
Bloch's theorem. I now briefly re-examine them from this 
point of view. Since the Hamiltonian H is periodic, an 
eigenbasis exists in which all eigenf unctions of H are of 
the form 

^n,k{r) = ^e*>„,fc(r), (2.6) 

with k G BZ and n a discrete index. The functions 
i'n,k{i') are ?7C-periodic and will be referred to as Bloch 
states. They satisfy H{k)ipn^k = ^n.k^n.k where e„^fc are 
the energies of the eigenstates. The effective Hamiltonian 
H{k) is given by 

H{k) = v(T-{-ihdr + hk + eA) + U = H+hva-k. (2.7) 

Normalizing 4'n,k{f) to unity over UC ensures that ^n,k 
is normalized to unity over the whole sample. The zero- 
eigenstates of Eq. (|2.5p are ?7C-periodic, and therefore 
correspond to fc = solutions in this labeling scheme. 
The corresponding normalized zero-energy eigenstates of 
H{k = 0) (with [/ = 0) are V±,o = f^*±,o- 

III. NEAR ZERO EIGENSTATES AND A 
NON-ZERO BUT SMALL SCALAR POTENTIAL 

The near-zero energy eigenstates in the vicinity of 
fe = can be studied by treating the term fiver -k in -ff (fc) 
perturbatively. The effect of a weak periodic scalar po- 
tential U can be treated perturbatively at the same time. 

To leading order in k and [/, the problem is solved by 
Bloch states of the form 

Urj,k{r) = X^,k'^+,o{r) + Xr^.k^-fiir), (3.1) 

where Xr,.fc = (x^,fc,X^,fc) satisfies h{k)x,,,k = £rj,kXn,k- 
Here rj = ± labels the bands (— for the valence band 
and + for the conduction band), e,,,fc is the energy of the 
state with label 77, k and h{k) is the Hamiltonian -ff (fc) 
projected onto the zero eigenspace. It has the form of a 
2x2 Dirac Hamiltonian in k-space: 

h{k)^n' ^ / dr^ il}^,Q{r)'^ [hv(T ■k + U{r)]il)^>^Q{r), 
Juc 

==> h{k) = hit a ■ k + fiaz (3.2) 



4 



where a term proportional to identity is omitted because 
it simply leads to a redefinition of the zero-energy. Here 
V — PN+N^v is the renormalized Fermi velocity, and P 
is the area of UC. The effective mass /i is given by 



/i 



UC 



(3.3) 



The signs of U{r) and F{r) are said to be correlated 
when there is a large overlap between regions of positive 
U and positive F and anti-correlated when there is a 
large overlap between regions of positive U and negative 
F. When the signs of U and F arc either correlated or 
anti-correlated, the two terms in the integrand tend not 
to cancel and is non-zero. 

The leading order perturbative expressions for the 
Bloch state u±,fc(r) in the vicinity of fe = is 



(3.4) 



with 6 = aig{kx + iky) and e±^k = ±-\/ {hv\k\)'^ + fi^ the 
energy of the state. 

It is worth emphasizing here that the same effective 
Hamiltonian h{k), and therefore the same mass term is 
induced in both valleys. This is to be contrasted with the 
mass term that arises when a staggered on-site potential 
is applied to otherwise clean graphene, by for instance 
placing the sample on an appropriate substrate.^ If the 
staggered potential takes on the value m on the A sub- 
lattice and — m on the B sublattice, then a mass term 
ma^ is induced in valley K while a mass term —maz is 
induced in valley K'^ 



IV. THE ANOMALOUS HALL EFFECT 

In the previous section I demonstrated that there is a 
gap in the spectrum and that a low energy description 
in terms massive Dirac fcrmions is possible. The gap has 
several important consequences. Firstly, low energy elec- 
trons can be localized in selected regions of a graphene 
sample by inducing a gap in regions where electrons are 
to be excluded. The gap can be controlled by electromag- 
netic fields that are smooth on the scale of the graphene 
lattice constant, and there is no need to break the sub- 
lattice inversion symmetryi** Secondly, since one has con- 
trol over the sign of the mass, it is possible to produce a 
sample in which the sign of the mass differs in different 
spatial regions. At the interface between two such re- 
gions, chiral edge states are expected to appear.^ Thirdly, 
as is the case in Ref. an electric field E — Ey applied 
to the bulk of the sample (the y-direction is chosen ar- 
bitrarily) will produce a quantized current jx = (JxyE 
in the direction perpendicular to it. This is remarkable 
since the magnetic field through the sample averages to 
zero. 

Here I briefly review the general theory for this 
phenomenon, which is known as the anomalous Hall 



effect ^i^J^iiiii^ii^iiiii^ It states that the Hall conduc- 
tivity axy may be non-zero, even in the absence of a 
net magnetic flux through the sample, provided time- 
reversal invariance is still broken. Furthermore it is 
quantized in units of e'^/h and can be written in terms 
of a set of integer topological invariants as axy = 
(e^/h) J2n Ch„^iiaiiia2ii3ii4j^ jjgj.g ^ jg ^j^g discrete in- 
dex that labels the energy bands and the sum ranges over 
all occupied bands. (The theory assumes that the Fermi 
energy is in a gap between bands so that there are no 
partially filled bands.) The invariant Ch„ is known as a 
Chern number, and can be calculated if the Bloch state 
for band n is known throughout the whole Brillouin zone. 
Since only the Bloch states in the vicinity of fc = and 
only for the two bands that touch at e = were found 
in Sec. IIIIl direct calculation of the Ch„ is not feasible 
here. 

Fortunately, there is a way to work around this prob- 
lem. The Chern numbers can only change when bands 
touchiii Before and after a touching of bands, they dif- 
fer by an integer. Only the Chern numbers of bands 
involved in the touching change. Furthermore, there is 
a sum-rule relating the Chern numbers before and after 
a band-touching. Suppose that, by varying a parame- 
ter in the Hamiltonian, one induces bands n and n + 1 
to touch before moving apart again. Let Ch^^ (Ch^) 
be the Chern numbers before (after) the bands touched. 
Then Ch^ + Ch<_^i = Ch> + Ch>_^i, i.e. the sum of the 
two Chern numbers is preserved. If bands n and n + 1 
are both occupied, axy stays constant although Ch„ and 
Ch„_|_i change. Thus axy only changes if the Fermi en- 
ergy lies in the gap between bands n and n+ 1. 

In order to calculate the change in Chern numbers, 
one only needs information about the bands that touch 
and only in regions of the Brillouin zone that are close to 
the points where the bands touch i^'^i^^ This information 
may be obtained by means of a perturbative expansion 
in the wave vector k and the external parameter that 
controls the closing of the gap. (The general procedure 
is the same as was used in Sec. IIIIl of this text.) If the 
first order coefficients do not vanish, the touching point 
is characterized by a Dirac fermion. Its effective 2x2 
Hamiltonian is brought into a standard form (essentially 
the valley-isotropic representation) by means of appro- 
priate 2x2 unitary transforms. It turns out that the 
change in Chern numbers is completely determined by 
the sign of the mass term in this representation. If the 
bands touch at several points in k space simultaneously, 
then each touching point gives an independent contribu- 
tion. 

For the graphene system I am considering, this leads 
to a result 



Ch> - Ch< = ±{sign[Ai>] - sign[//<]}. 



(4.1) 



Here + (— ) refers to the conduction (valence) band and 
fi^ (/^^) is the mass induced by the scalar potential be- 
fore (after) the bands touched. There are two contribu- 
tions to the result, one from each touching point. Since 
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the same mass term appears at both touching points 
(in the valley- isotropic representation) , the contributions 
from the two touching points are the same. If the sign 
of the mass were opposite in the two valleys, as is the 
case for a substrate-induced mass, the two contributions 
would have canceled. 

There are two ways to change the sign of the mass 
term [cf. Eq. (|3.3p ]. One may either change the sign 
of the external potential ?7, or one may change the sign 
of the magnetic field and hence of F. According to the 
theory, both methods give the same change in the Hall 
conductivity, Aa^y = Aa^yiU —>■ —U) = Aaxy{B — > 
—B). The transformation B — > ~B is equivalent to time- 
reversal, which changes the sign of the Hall conductivity: 
o'xy = —<7xy- Therefore AcFxy = 'ioxy- Thus, knowing 
the change in <Jxy when the mass changes sign, one also 
knows the actual value of Oxy 

If is the Hall conductivity of a sample in which 
the mass term has sign ± then a^y — Oxy = ~2e'^ /h. 
Therefore, according to the above argument 

V. ORIGIN OF NON-ZERO HALL 
CONDUCTIVITY 

Suppose that the topological theory discussed in the 
previous section was not known. Would it have been 
possible to show that the graphene system studied here 
displays a Hall effect, without at the same time essen- 
tially deriving the topological theory? As I show in this 
section, the answer is 'yes'. In the process, the underlying 
physical mechanism that here produces the anomalous 
Hall effect is identified. 

The starting point of the argument is to consider the 
two zero-energy solutions 5'±,o e^^ that are obtained 
when U — 0. A key feature of these wave functions 
is that the sign of F is anti-correlated with that of 
_B i.e. is negative on average where the magnetic 
field points in the positive z direction and positive on 
average where the magnetic field points in the nega- 
tive z direction. To prove this assertion I show that 
the number P = (e/h) J^^ dr^ B{r)F{r) is negative. 
The proof is elementary. From Eq. (|2.3p follows that 
P = Jjj^ dr'^ B{r)d^B{r). Using integration by parts, 
and invoking Gauss's theorem I find 

P=(f) drj_ BdrB- dr^ldrBl"^ <0. (5.1) 

JdUC JuC 

^ V ' 

=0 

The integration in the first term is around the border 
of UC and dr± is normal to the border (in the x ~ y 
plane) pointing outwards and yields zero because B is 
C/C-periodic. The result [Eq. (|5.ip ] implies that a parti- 
cle in the state \E'_,o (^+,o) is largely confined to regions 
where B points in the positive (negative) z direction. In 



Sec. HI] it was noted that the eigenstates ^'±,o, have op- 
posite sublattice polarization. The correlation between 
B and F then implies that sublattice polarization is cor- 
related with the sign of the magnetic field. For instance, 
in valley if, regions where the magnetic field points in 
the negative z direction are A sublattice polarized while 
regions where the magnetic field points in the positive z 
direction are B sublattice polarized. This is consistent 
with the numerical results presented in Figure 3 of Ref. 
[TtI . There it was pointed out that the state VP-^o closely 
resembles a state in the zero Landau level of a graphene 
sample with a constant magnetic field that points in the 
positive z direction. This would have explained the exis- 
tence of a Hall effect, were it not for the state ^'+,o that 
also belongs to the zero eigenspace. This state is similar 
to a zero Landau level state in a magnetic field pointing 
in the negative z direction, and it is responsible for a Hall 
effect that cancels that of 5'-,o- 

In order to have a Hall effect a perturbation that lifts 
the degeneracy between 5'+^o and 5'_.o is needed. This is 
the role of the electrostatic potential U . If the sign of U 
is anti-correlated with that of B it will raise the energy 
of the state 5'+,o while lowering the energy of the state 
\l/_,o. This leads to a gap with associated positive mass. 
A scalar potential U whose spatial profile is correlated 
with that of B will do the opposite and leads to gap 
associated with a negative mass. 

Turning to the expression [Eq. (|3.4] for the U ^ 
valence band Bloch states one sees that the dom- 

inant term in the linear combination is proportional to 
g-sign(Ai)-F result, valence band states are confined 

to regions of positive B when U is anti-correlated with 
_B, and confined to regions of negative B when U is cor- 
related with B. This is consistent with the response of 
^1^-1- to U discussed above. 

One thus arrives at the fundamental result that there 
is a finite positive (negative) magnetic flux through the 
region of space occupied valence band electrons when the 
mass is positive (negative). The valence band will then 
experience a 'usual' Hall effect. An electric field E = Ey 
will, through the Lorentz force, lead to a current in the 
direction — sign(/i) x. Electrons in the conduction band 
would have seen an opposite magnetic field and produced 
an opposite Hall effect. However, the Fermi energy is in 
the gap between the two bands so that no conduction 
band states are occupied. The sign obtained here for 
the Hall current agrees with the result derived from the 
topological theory [cf . Eq. (|4.2p ] . 

So far I have said nothing about the quantization of the 
Hall current. Indeed, the above argument does not prove 
that the anomalous Hall conductivity precisely equals 
=pe^//i. It only indicates that Uxy is finite and predicts 
its sign. However, the argument does relate the anoma- 
lous Hall effect in graphene with periodic electromagnetic 
fields to the usual integer Hall effect in a constant mag- 
netic field. It thus shows that the quantization of axy 
has the same origin in both cases. 

Finally, a word on the nanotube system: The Hall ef- 
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feet found there^^ is not produeed by states loealized 
where the magnetic field has a well-defined direction, al- 
though such states are of course also present in that sys- 
tem. Rather it is produced by higher energy states that 
are localized where the magnetic changes sign. Thus, a 
different mechanism is responsible for the effect reported 
in Ref. [H. Other features also distinguish it from the 
Hall effect considered in the present work: ft is not ac- 
companied by a gap between the conduction and valence 
bands. It is associated with two sets of states on oppo- 
site sides of the circumference of the nanotube that carry 
current along the length of the nanotube but in opposite 
directions. The quantization is approximate. This is to 
be contrasted with the following properties of Hall effect 
discussed in the present work. It appears when the Fermi 
energy is in the gap between the conduction and valence 
bands. All the current travels in one direction through- 
out the whole sample. The Hall quantization of the Hall 
conductivity is sharp. 

The conclusion is that the Hall effect discussed in this 
paper has not yet been found in the context of carbon 
nanotubes. But is it there? The topological theory can- 
not answer this question. It only works for systems that 
are much larger than the periodicity of the magnetic field, 
whereas the length of the nanotube system in the circum- 
ferential direction is the same as the periodicity of the 
magnetic field. It is not clear if the intuitive argument 
presented in this section can be applied either. It also in- 
volves viewing the system with a spatial resolution that 
is larger than the periodicity of the magnetic field. The 
question therefore remains open. 



VI. EXAMPLES 

Example 1: As a simple example that illustrates the 
general results derived in this text, consider the magnetic 
field B{r) = Bq (sin27ra;/A -I- sin27ry/X) with periodicity 
A in both the x and y directions. 

It is convenient to define the magnetic length 1^ — 
y^h/e\BQ\ and a dimensionless constant (3 = eX^lBol/h = 
{.X/lmY- According to Eq. (|2.4p . F is given by 
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(6.1) 



This gives the clearest possible illustration that 5'±,o are 
localized to regions of positive or negative B. The degree 
of localization is controlled by the ratio X/lm- 

The renormalized Fermi velocity can be calculated ex- 
plicitly as 

(6.2) 



/o(/3/47r2)2 



where /q is the zero'th modified Bessel function of the 
first kind. When /? ;2> 1 or equivalently /,„ <C A, the 
asymptotic form v ~ (3e~^^'^^ v jl'K is valid. In this limit 
the Fermi velocity becomes very small and bands become 



nearly dispersionless in the region of fc = 0. The result is 
very similar to the nanotube result 3^ This is due to the 
fact that the spatial profile of the magnetic field in the 
nanotube case is also sinusoidal. In Figure [1] the solid 
line shows the renormalized Fermi velocity (solid line) 
plotted as a function of the magnetic field strength Bq. 

As was pointed out in Ref. lf| and [TtI the flat bands 
are reminiscent of Landau levels. It should however be 
pointed out that flat bands do not guarantee a Hall effect: 
The non-relativistic Schrodinger equation in the presence 
of the magnetic field of this example was studied numer- 
ically in Ref. Flat bands were found in the limit 
A/^m ^ 1- For that system the Chern numbers of flat 
bands were found to be zero, leading to a zero Hall con- 
ductivity. 

Neither are the flat bands a necessary condition for 
the Hall effect to occur. The effect also occurs for 1^ > 
A when the renormalized Fermi velocity is of the same 
order of magnitude as the unrenormalized Fermi velocity. 
Counter-intuitively, it is not the flat bands, but the Dirac 
cones at the touching points between bands, that lead to 
the Hall effect. 
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FIG. 1: The renormalized Fermi velocity v (solid line) ac- 
cording to Eq. (|6.2|l and the mass ^ (dashed line) accord- 
ing to Eq. (|6.3|l for the model with magnetic field B{r) = 
Bo {sm2Tix / \ + sin2Tiy / \) perpendicular to the graphene 
plane. 



Returning to the analysis of the example I now intro- 
duce a scalar potential U{r) = —UoB{r)/BQ with the 
same spatial profile as the magnetic field. The induced 



/i = sign(5o)[/o 



/i(/3/4^2 
/o(/3/4^2 



sign(Bo)t/o, (6.3) 



where Ii is the first modified Bessel function of the first 
kind. The asymptotic form is valid in the same limit 
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as before, namely when Im X. This result is straight 
forward to interpret. In the strong magnetic field hmit 
the eigenstates become exponentially well-confined to the 
maxima and minima of the magnetic field. These regions 
are also the maxima and minima oi U , so that the effec- 
tive potential seen by the eigenstates has magnitude \ Uo\. 
The behavior of the mass as a function of the magnetic 
field strength Bq is shown (dashed line) in Figure [TJ 

Example 2: Finally, I consider an example in which the 
magnetic and electric fields are constant in one direction 
and periodic in the other. The full spectrum can be ob- 
tained analytically. The approximate formulas for the 
low energy spectrum in the previous sections can there- 
fore be tested against exact results. 

The magnetic field and electrostatic potential are taken 

as 



nez 

u - -UoX J2 [^(y - - ^(y - + 1)^)] •(6-4) 



nez 



The function F is given by 
A2 /I 



F = 



211 



y mod A 



A 



(6.5) 



According to the results presented in Sec. IIIIl the renor- 
malized Fermi velocity and mass are 



4sinh(/3/4) 



sign(Bo) — , 



(6.6) 



with (3 defined as in Example 1. Again the renormalized 
Fermi velocity is an exponentially decreasing function of 
P and hence of |i?o|- 

The full spectrum of the model can be obtained from 
the transcendental equation 



^ 2 ' 

sin(^^) sin(^l^) 



cos(fcyA) = cos(^^) cos(^|— ) 



P+ 



P- 



(6,7) 



where p± — ^ (e/hv)^ — fc±, k± — ± A/2Z,^j, e is the 
energy and the wave vector is fc = k^x + kyij. In Fig- 
ure the exact dispersion of the two bands closest to 
e = is plotted in the x and y directions. On the 
same plot is also shown the hyperbolic dispersion rela- 
tion e = ±-\/ (?iu|fc|)^ -I- fj?^ with mu and v as calculated 
in Eq. (|6.6p . It is seen that the low energy dispersion 
relation is reproduced very well by the approximate for- 
mulas. 



VII. CONCLUSION 

In this paper I showed analytically that when a 
jraphene sample is exposed to a periodic magnetic field 
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FIG. 2: The low energy dispersion of the Dirac equation with 
the magnetic field and electrostatic potential of Eq. (|6.4|l . 
The dashed line shows the exact result for the wave vector k 
chosen along the x direction. The dot-dashed line shows the 
exact result for the wave vector fc along the y direction. The 
solid line shows the Dirac dispersion e = {fiv\k\)'^ -\- 
with /i and v as in Eq. (|6.6p . A magnetic length Im = A/2 
and a potential strength Uo — hv/AX was used. 



that (a) is smooth on the scale of the graphene lattice, 
and (b) is zero on average, the valence and conduction 
bands still touch at two Dirac points, but the Fermi ve- 
locity is reduced. The reduction the same in all direc- 
tions and is exponential in the magnetic field strength. 
For simple models that introduce no other length scales 
beyond the magnetic length 1^ and the periodicity A of 
the magnetic field, the renormalized Fermi velocity v be- 
haves as v/v oc exp(cA^/Z^) with c of order unity. The 
zero eigenspace is two-fold degenerate. A basis exists in 
which one eigenstate confines particles to regions where 
magnetic field points in the positive z direction and to 
sublattice A (B) in valley K, (K'). The other eigenstate 
does the opposite, i.e. confines particles to regions where 
the magnetic field points in the negative z direction and 
to sublattice B (A) in valley K, {K'). An electrostatic 
potential that is either correlated or anti-correlated with 
the magnetic field induces a gap between the valence and 
conduction bands. These analytical results are consistent 
with previously obtained numerical results.'^ I showed 
that the gapped state supports a quantized Hall effect 
with axy = ±e^//i. The positive sign refers to the cor- 
related case and the — sign refers to the anti-correlated 
case. This is an instance of the known phenomenon of 
a Hall effect in zero overall flux.— It can be explained in 
terms of the topological properties of Bloch states of the 
occupied bands4'^ i^°'^^i^^'^'^i^^'^^ In the present system I 
showed that the effect is also simply related to the usual 
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integer Hall effect in a non-zero average magnetic field: Acknowledgments 
In the presence of the electrostatic potential, the spatial 
region in which the filled valence band states are local- 
ized, is permeated by a finite total magnetic flux. As 
a result, the valence band electrons see a magnetic field 

with a non-zero average even though the magnetic field This research was supported by the National Research 
averages to zero over the sample as a whole. Foundation (NRF) of South Africa. 



* Electronic address: "isnyman@sun.ac.za' 

^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 

^ M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 

Phys. 2, 620 (2006). 
^ C. W. J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008). 
^ S. Y. Zhou, G.-H. Gweon, A. V. Fedorov, P. N. First, W. 

A. de Heer, D.-H. Lee, F. Guinea, A. H. Castro Neto, and 

A. Lanzara, Nature Mat. 6, 770, (2007). 

^ F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
® F. D. M. Haldane, Phys. Rev. Lett. 93, 206602 (2004). 

C. L Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005). 

® R. Karplus and J. M. Luttinger, Phys. Rev. 95, 1154 
(1954). 

^ D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. 

den Nijs, Phys. Rev. Lett. 49, 405 (1982) 
^° Q. Niu, D. J. Thouless, and Y. Wu, Phys. Rev. B 31, 3372 
(1985). 

J. E. Avron, R. Seller, and B. Simon, Phys. Rev. Lett. 51, 
51 (1983). 

B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 

" M. Kohmoto, Ann. Phys. (N.Y.) 160, 343 (1984). 
" M. Oshikawa, Phys. Rev. B 50, 17357 (1994). 

M. Onoda and N. Nagaosa, J. Phys. Soc. Jpn. 71, 19 

(2002). 



F. Guinea, M. I. Katsnelson and M. A. H. Vozmediano, 
Phys. Rev. B 77, 075422 (2008). 
^"^ T. O. Wehling, A. V. Balatsky, A. M. Tsvelik, M. L Kat- 
snelson, and A. I. Lichtenstein, Euro. Phys. Lett. 84, 17003 
(2008). 

^® A. Krakovsky Phys. Rev. B 53, 8469 (1996). 

M. Taillefumier, V. K. Dugaev, B. Canals, C. Lacroix, and 
P. Bruno, Phys. Rev. B 78, 155330 (2008). 
Unfortunately 'anomalous Hall effect' is also the standard 
terminology for another phenomenon in graphene, namely 
the unusual quantization of Landau levels. This is not what 
the term refers to in this text. Rather it refers to a Hall 
effect in a zero average magnetic field. 
L. Brey and H. A. Fertig, arXiv:0904.0540 (2009). 

^2 H.-W. Lee and D. S. Novikov, Phys. Rev. B 68, 155402 
(2003). 

S. Bellucci, J. Gonzalez, F. Guinea, P. Onorato, and E. 
Perfetto, J. Phys.: Condens. Matter 19, 395017 (2007). 
R. Jackiw, Phys. Rev. Lett. 29, 2375 (1984). 
2^ J. K. Pachos and M. Stone, Int. J. Mod. Phys. B 21, 5113, 
(2007). 

^® The sign of the mass term in a given valley depends on the 
representation. All statements made in this text hold for 
the valley-isotropic representation. Care should be taken 
when comparing with results in the literature since several 
other representations are also common. See footnote 2 of 
Ref.i. 



